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Abstract 

We show that a hierarchical CPT violating neutrino spectrum can simultane- 
ously accommodate all published neutrino data, including the LSND result, with a 
total x 2 which is almost identical to the CPT conserving best fit. In our scenario 
the oscillation signal seen by the KamLAND experiment in antineutrinos is inde- 
pendent of the LMA solar oscillation signal seen in neutrinos. A larger antineutrino 
mass splitting accounts for the LSND signal and also contributes to atmospheric 
oscillations. Because of this feature, there is some tension between our model and 
certain Super-K atmospheric results. Thus, if LSND did not exist, our model would 
survive only at the 99% confidence level; with LSND, our model is (essentially) 
statistically equivalent to the CPT conserving best fit. 



1 Introduction 



CPT violating neutrino masses allow the possibility PQ - [3] of reconciling the LSND [5], 
atmospheric jB], and solar oscillation j7[|Hj data without resorting to sterile neutrinos. As 
argued in [2], there are good reasons to imagine that CPT violating dynamics couples 
directly to the neutrino sector, but not to other Standard Model degrees of freedom. An 
explicit CPT violating model of this type was presented in [I]. 

KamLAND j^j, a medium baseline reactor antineutrino disappearance experiment, is sen- 
sitive to antineutrino mass-squared splittings in the 10~ 4 eV 2 range characteristic of the 
large mixing angle (LMA) solar neutrino scenario. The KamLAND collaboration has 
recently reported ^0] an electron antineutrino survival probability which is significantly 
less than one: 

P(V e -> V e ) = 0.611 ± 0.085 ± 0.041 . (1.1) 

If the neutrino mass spectrum conserves CPT, then this result is consistent with the 
LMA interpretation of solar neutrino oscillations. If the neutrino mass spectrum violates 
CPT, however, the KamLAND result provides no information about solar oscillations, 
but rather constrains the splittings in the antineutrino spectrum. 

In this paper we show that a hierarchical CPT violating neutrino spectrum can simultane- 
ously accommodate the oscillation data from LSND, atmospheric, solar and KamLAND, 
as well as the nonobservation of antineutrino disappearance in short baseline reactor 
experiments. Comparing our model to the global neutrino data set we obtain a total 
X 2 = 201 for 228 degrees of freedom. This can be compared to the usual CPT conserving 
best fit, which has a total x 2 = 199 for 232 degrees of freedom [TTj . 

While the total x 2 are about the same, there are striking differences when we make a 
more detailed comparison to the data. The CPT conserving best fit parameter set (two 
mass differences and 3 mixing angles) matches remarkably well to the global data set, 
with the glaring exception of the LSND result. This single data point contributes a A% 2 
of about +12, leading to widespread speculation that the LSND result is simply wrong, 
and will be disconfirmed by MiniBooNE. 
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As we will see, the hypothesis of a CPT violating neutrino sector (four mass differences 
and 6 mixing angles) leads to a completely different conclusion. In this framework our 
(non-optimized) parameter choices give an excellent fit to LSND. The largest single dis- 
crepancy with the global data set instead occurs with the atmospheric Super-K multi-GeV 
muons and thru-going muons. In either case the A% 2 is no worse than +4. Thus instead 
of casting a skeptical eye upon LSND one is led to speculate that the forthcoming re- 
analysis of the Super-K atmospheric data will improve this discrepancy. Indeed the new 
fluxes announced for the Super-K re-analysis definitely work in this direction. 

We also make a dramatic prediction for the observation of atmospheric neutrinos using 
the MINOS detector [T3]. Because the MINOS detector discriminates positive and neg- 
ative charge, this experiment can disentangle the neutrino and antineutrino components 
of atmospheric oscillations in a straightforward way. As the mass differences in the atmo- 
spheric sectors differ by orders of magnitude in our scenario, MINOS will be able to tell 
them apart easily. 

2 The spectrum 

To analyze all the possible CPT violating spectra is not an easy job. With four mass dif- 
ferences and six mixing angles (not taking into account the two CP violating phases which 
participate in oscillations) a complete scan of the whole parameter space is impractical. 
However, thanks to the available experimental data, it is possible to reduce the allowed 
regions to two sets of well-differentiated spectra with (quasi) orthogonal experimental 
signatures. 

The easiest way to make contact with the experimental results is in terms of the neutrino 
survival and transition probabilities, which are given by 



P(v a -> vp) = 5 af3 - 4 UaiUpiUajUpj sin 
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for neutrinos and 

— - - TAmf.L" , 

(2.2) 



P(y a -> v p ) = 5 aP -A J2 U ai U piU aj U Pj sin 5 

i>j=i 
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for antineutrinos. Here the matrix U = {U ai } (U = jfJoi}) describes the weak inter- 
action neutrino (antineutrino) states, v a , in terms of the neutrino (antineutrino) mass 
eigenstates, j^. That is, 

v a = Yl U ™ v i and v <* = ^ZU ai Vi (2.3) 

i i 

where we have ignored the possible CP violation phases in both matrices and took them 
to be real. The matrices can be parametrized as follows: 

' C12C13 S12C13 S13 \ 

U = -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 (2.4) 

V S12S23 — C12C23S13 — C12S23 — S12C23S13 C23C13 j 

and similarly for U. In Eq. (12. 1|) L denotes the neutrino flight path, i.e. the distance 
between the neutrino source and the detector, and E is the energy of the neutrino in the 
laboratory system. 

Regarding the mass spectrum of the three neutrinos we assume that it is hierarchical and 
thus characterized by two different squared masses 

a2 2 2 JA2 2 2 

Am 12 = m 2 — m 1 and Am 13 = m 3 — m 1 

whose numerical values are rather different, i.e. Am^ ^> Am^ 2 and similarly for the 
antineutrinos. Having said that, it becomes apparent that the larger mass-squared differ- 
ence in the neutrino sector will be related to the atmospheric neutrino signal observed by 
Super-Kamiokande, while the smaller one will drive the solar neutrino oscillations. In the 
antineutrino sector, the largest mass difference will provide an explanation to the signal 
observed in LSND, while the smaller one is the one which might have been (mis) identified 
by KamLAND ClS db confirmation of LMA. 

The key ingredient to sort out the antineutrino spectra are reactor experiments. Their 
results indicate fSl ^ na ^ electron antineutrinos produced in reactors remain electron 
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Figure 1: Possible neutrino mass spectrum with almost all the electron content in the heavy 
state. Although the figure shows an example of large mixing, our approach is agnostic about 
the mixing matrix. The flavor content is distributed as follows: electron flavor (red), muon 
flavor (brown) and tau flavor (yellow) 



antineutrinos on short baselines. As the distance traveled by our antineutrinos is small we 
can forget about the smallest mass difference and average the other two, thus the survival 
probability can be expressed as 

P(V e -+17 e ) = l-2U 2 e3 (l-Ul 3 ) . (2.5) 

It is clear that there are two possible ways to achieve a survival probability close to one, 
i.e. C/ e 3 can be almost one or almost zero. Physically this means that we can choose 
between having almost all the antielectron flavor in the heavy state (or in the furthest 
away state) or just leave in this state almost no antielectron flavor. The first possibility 
(which is depicted in Fig. 1) is the one we explored in our previous works. This spectrum 
predicts for KamLAND a survival probability consistent with one. Since this is strongly 
disfavored by the KamLAND result , we instead pursue the second possibility, which 
is represented by the spectrum shown in Figure 2. 

This second family of spectra is characterized by a strong violation of CPT in the mass 
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Figure 2: Possible neutrino mass spectrum with almost no electron content in the heavy 
state. Although the figure shows an explict mixing pattern, there is a whole family of 
mixing matrices that can do an equally good job. The flavor content is distributed as 
follows: electron flavor (red), muon flavor (brown) and tau flavor (yellow) 

differences but a much slighter effect in the mixing matrix. This is seen in Fig. 2 where 
the flavor distribution in the neutrino and antineutrino spectra is rather similar. The 
most distinctive feature of this family of solutions is its #23, which lives far away from 
maximal mixing, or in other words which has a large component of antitau neutrino in the 
heavy state. The small antimuon neutrino component in the heavy state is not bounded 
by the non observation of muon neutrino disappearance over short baselines in the CDHS 
experiment as the antineutrino component in this experiment was minimal. 

KamLAND could have observed an oscillation signal driven by the smaller antineutrino 
mass splitting and interpreted it as LMA oscillations. To explicitly see how this might 
have happened, we will choose two sample points in our parameter space and calculate the 
transition probabilities for it. Let us emphasize that we have not performed a chi-squared 
fit and therefore the points we are selecting (by eye and not by chi) are not optimized to 
give the best fit to the existing data. Instead, they must be regarded as two among the 
many equally good sons in this family of solutions. 
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The point we have chosen has 6 13 = .08 , 6 23 = .5 , 9 12 = .6 , Am 2 2 = 5 • 1CT 4 eV 2 
and Amjj = Oil) eV 2 . Since we are dealing with an antineutrino signal, we do not need 
to identify either the flavor distribution or the mass eigenstates of the neutrino sector. 
We will do it later, when showing the zenith angle dependence this model predicts for 
Super-Kamiokande atmospheric neutrinos. 



The survival probability measured by KamLAND is given by 

"Am 2 3 L 



KamLAND 



= 1-4XTJ1-U 



e3 



sin 
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4U 2 el U 2 e2 sm 2 



Am 2 2 L 
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(2.6) 



where the second term (proportional to U e3 ) is negligible. Plugging our numbers in, it is 
straightforward to see that Poland ~ -6 regardless of whether the mass difference that 
drives the solar neutrino oscillations belongs to the LMA region. 

By the same token, we can calculate the probability associated with the LSND signal. It 
is given by 



-Plsnd = 4^ 3 ^ 3 sin 2 



Amj 3 L 
AE 



(2.7) 



where we have neglected terms proportional to Am 2 2 which are irrelevant for such small 
distances. As the reader can easily verify, we predict a P LSND ~ .0022 in excellent agree- 
ment with the LSND final analysis: 



LSND-final 



0.00264 ± .00081 . 



(2.8) 



The only piece of experimental evidence involving antineutrinos which remains to be 
checked is the signal found for Super-Kamiokande atmospheric neutrinos. As we are 
introducing an antineutrino mass difference roughly two orders of magnitude larger than 
the Super-K best fit point (for an analysis with two generations and conserving CPT), 
there is cause for concern. In fact we pass this test as successfully as we did the others. To 
see this, we have first to state the parameters in the neutrino sector. Once more they have 
been chosen almost randomly from the different analyses available in the literature and 
are given by 13 = .08 , 6 23 = .78 , 6 12 = .52 , Amj 2 = 1 • 10~ 4 eV 2 and Am 2 13 = 2.8 • 10~ 3 
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eV 2 . We stress that although we have chosen a point in the LMA region, the particular 
election of both Amf 2 and 612 does not affect the quality of the agreement with the data. 

With these parameters we have calculated the zenith angle dependence of the ratio (ob- 
served/expected in the no oscillation case) for muon and electron neutrinos for the sub- 
GeV and multi-GeV energy ranges (remember that since Super-K is a water Cherenkov 
detector it does not distinguish neutrinos from antineutrinos and washes out any possible 
difference between the conjugated channels). The results are shown in Fig. 3 where we 
have also included the experimental data for the sake of comparison. As we have closely 
followed the spirit of the calculation in fOllH], we refer the reader to this article for de- 
tails and skip the technicalities. We worked in a complete three generation framework 
and included matter effects. 

In Fig. 4 we show the comparison to Super-K for our second example point. For this point 
we have chosen 6 13 = .08 , 6 23 = .5 , 6 12 = .785 , Am 2 2 = 7 • KT 5 eV 2 and Am 2 13 = 0(1) 
eV 2 . Note that this point is consistent with the best-fit point of KamLAND [TU] . 

3 Comparison with other analyses 

It is clear from Figs. 3 and 4 that our CPT violating spectrum does pretty well, with a 
moderate discrepancy apparent only for the multi-GeV muons. In order to understand 
the results it is important to remember that due to production and cross section effects 
Super-K is dominated by neutrinos, with antineutrinos a minor (but not negligible) con- 
tribution. This is only part of the story, however, since an analysis done by the Super-K 
collaboration allowing for CPT violation did not allow (at 99% C.L.) a mass difference in 
the antineutrino sector as drastically different as the one we are proposing. The difference 
is due to the fact that the Super-K analysis was done in a two generation context and, 
most importantly, forcing the two mixing angles to be maximal. This latter fact indeed 
maximizes the antineutrino contribution and compels the antineutrino mass difference to 
take the closest possible value to the neutrino one. By the same token, if we want a large 
antineutrino mass difference, we expect to improve the fit to Super-K data by combining 
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this with a non-maximal mixing angle, which suppresses the antineutrino contribution to 
the Super-K oscillation signal. 

The two vs three generation analysis has also an impact, as is seen by inspecting the 
transition probability for muon antineutrinos into tau antineutrinos, which is given by, 

From this formula it becomes apparent that for neutrinos coming from above only the 
largest mass difference contributes. However, for those neutrinos which have travelled 
through sizeable portions of the Earth and have covered distances of the order of 10 4 km, 
the second mass difference also plays a role. This contribution (which does affect the final 
result, especially for sub-GeV neutrinos) is neglected if only one mass difference is taken 
into account. 

Our analysis agrees with the spirit of the findings in Ref where a two generation 
approximation that didn't include matter effects was used. Also a simplified analysis based 
only on the up/down asymmetry in the number of multi-GeV events (in the CPT violating 
case) is available in the literature |18j . which used an older Super-K data set. If one uses 
(as we do) the result from the full 1490 day of SK-I data, i.e. = -.288 ± .030 [H| the 
CPT violating case (which gives for the sample points we have being using = —.27) 
turns out to be as acceptable as the CPT conserving one (with A^ = —.33 for maximal 
mixing but where smaller values can be obtained by departures of maximal mixing). In 
all the cases the electron neutrino asymmetry is consistent (within experimental errors) 
with zero. 

In order to get a quantitative estimate of how well our CPT violating spectrum fits the 
global data set not including LSND, we have compared (for both of our sample points) 
the total x 2 to the minimum total \ 2 which is obtained for different values of the 10 
input parameters. Fixing #13 = #13 = .08 our (admittedly crude) program finds that the 
minimum total x 2 occures for #23 = #23 = -8, 6>i 2 = 12 = .6, Am^ = Am^ = 6 ■ 10~ 5 
eV 2 , and Am 2 3 = Am 2 3 = 3 • 10~ 3 eV 2 . This agrees rather well with more sophisticated 
best fit analyses. We find the A% 2 for either of our sample points to be smaller than 9.6. 



Am 2 3 L 
AE 



- AU^U^U^U^sm 



Am| 2 L 
AE 



(3.1) 
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Making the conservative observation that variations in 13 , 13 don't affect the total x 2 
as long as the central values are small, we can interpret A\ 2 > 9.6 for 8 parameters as 
occurring with probability 29% due to Gaussian fluctuations in the global data set. 

However this result is questionable since, as it happens, the minimum total x 2 is obtained 
for values of the 8 input parameters which are very close to CPT conserving. Thus one 
can argue that we are only interested in contributions to A% 2 arising from fluctuations 
in the data which simulate variation in the 4 independent CPT violating parameters. In 
that case we are dealing with A% 2 > 9.6 for 4 parameters, which occurs with probability 
5%. 

After our model was proposed it was pointed out by Gonzalez-Garcia, Maltoni and 
Schwetz [TI] that the comparison with Super-K atmospheric data needs to include the 
zenith angle dependence of the thru-going muons. Because the data contains too many 
events near the horizon compared to the prediction of our model, this turns out to give 
a large contribution to the A% 2 in an analysis like that described above. Using a CPT 
violating spectrum very close to our sample points, these authors obtain a A% 2 = 12.7, 
which for 4 parameters implies a fluctuation with probability 1.3%. This is not especially 
encouraging, although it is worth pointing out that the Standard Model fit to the global 
electroweak data set has a confidence level of only CL = 0.02 [2*U] . 

Thus, if the LSND experiment did not exist, we would have concluded that our model 
has a confidence level of CL = 0.01. The fact that it competes well with the CPT 
conserving scenario once LSND is included is an indication that LSND is not statistically 
very compatible with the rest of the global data set. Indeed Gonzalez-Garcia, Maltoni 
and Schwetz have estimated, using their own customized goodness-of-fit test, that the 
LSND and all-but-LSND data sets are only statistically compatible at confidence level 
CL = 7.5 x 10~ 4 . It is important to note that this stringent result, although it employs 
the CPT violating parameter space, does not represent the confidence level for our model. 
As seen above, the confidence level for our model if LSND did not exist would be CL = 
0.01. With the entire global data set including LSND, the confidence level is CL = 0.74 
(A^ 2 > 2 for 4 parameters). 
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4 Outlook 



Once we have established that a CPT violating mass spectrum as the one shown in Fig. 
2 can account for all the available experimental evidence data, it is time to ask how we 
might confirm CPT violation in future data. 

The most straightforward answer is through experiments able to run in both modes (neu- 
trino and antineutrino), by simple comparison of the conjugated channels. The first of 
them is MiniBooNE, which is meant to close the discussion about LSND one way or the 
other. MiniBooNE is taking data and is expected to give a definite answer to the CPT 
question after some years of running. Needless to say we expect MiniBooNE to confirm 
LSND only when running in the antineutrino mode. 

For our type of spectrum, the observation of atmospheric neutrinos using the MINOS 
detector [T2j is also ideal. Because the MINOS detector discriminates positive and neg- 
ative charge, this experiment can disentangle the neutrino and antineutrino components 
of atmospheric oscillations in a straightforward way. As the mass differences in the atmo- 
spheric sectors differ by orders of magnitude in our scenario, MINOS will be able to tell 
them apart easily. 

A positive oscillation signal at KamLAND (which could be a misidentification of a CPT 
violating spectrum as LMA) and Borexino [21] finding a day /night asymmetry (evidence 
of a LOW solution [22]) or a seasonal variation (an indication of VAC (221) wm P om t 
towards CPT violation. Indeed a conflict between KamLAND and Borexino results would 
constitute strong evidence for CPT violation even if LSND is discontinued by MiniBooNE. 
It is worth mentioning that, despite the accumulation of strong evidence for the LMA 
solution, some interesting discrepancies remain j2S]- 

All in all, CPT violation has the potential to explain all the existing evidence about 
neutrinos with oscillations to active flavors. Such a scenario makes distinctive predictions 
that will be tested in the present round of neutrino experiments. One should always bear 
in mind that so far we have no evidence of CPT conservation in the neutrino sector. The 
true status of CPT in the neutrino sector might be established by the combined results of 
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KamLAND, Borexino and SNO, and certainly by MiniBooNE. In the atmospheric sector 
MINOS is the ideal experiment for such a test. 



Acknowledgments 

We are grateful to Andre de Gouvea, Bill Louis, Michele Maltoni, and Steve Mrenna for 
comments and assistance. This research was supported by the U.S. Department of Energy 
Grant DE-AC02-76CHO3000. 



References 

[1] H. Murayama and T. Yanagida, Phys. Lett. B 520, 263 (2001). 

[2] G. Barenboim, L. Borissov, J. Lykken and A. Y. Smirnov, JHEP 0210, 001 (2002) 
|arXiv:hep-ph/0108199| . 

[3] G. Barenboim, L. Borissov and J. Lykken, Phys. Lett. B 534, 106 (2002) 
|arXiv:hep-ph/02010"80| . 

[4] G. Barenboim and J. Lykken, Phys. Lett. B 554, 73 (2003) |arXiv:hep-ph/0210411 . 

[5] LSND Coll. (C. Athanassopoulos et al), Phys. Rev. Lett. 75, 2650 (1995); 77, 3082 
(1996); 81, 1774 (1998); W.C. Louis for the LSND Coll., Nucl. Phys. B (Proc. Suppl.) 
91, 198 (2001); LSND Coll. (A. Aguilar et al), |arXiv:hep-ex/0104049l 

[6] Super-Kamiokande Coll. (Y. Fukuda et al.), Phys. Lett. B433, 9 (1998); Phys. Rev. 
Lett. 81,1562 (1998); 82, 2644 (1999); Phys. Lett. B467, 185 (1999). 

[7] Homestake Coll. (B.T. Cleveland et al), Astrophys. J. 496, 505 (1998); GALLEX 
Coll. (W. Hampel et al), Phys. Lett. B447, 127 (1999); V.N. Gavrin, for the SAGE 
Coll., Nucl. Phys. B (Proc. Suppl.) 91, 36 (2001); E. Bellotti, for the GNO Coll., 
Nucl. Phys. B (Proc. Suppl.) 91, 44 (2001); SuperKamiokande Coll. (S. Fukuda, et 



12 



al.) Phys. Rev. Lett. 86, 5651 (2001); 86 5656 (2001); SNO Coll. (Q.R. Ahmad et 
al), Phys. Rev. Lett. 87, 071301 (2001); 89, 011301 (2002). 

[8] SNO Coll. (Q.R. Ahmad et al), Phys. Rev. Lett. 89, 011302 (2002). 

[9] A. Piepke [KamLAND Collaboration], Nucl. Phys. Proc. Suppl. 91, 99 (2001); 
S. A. Dazeley [KamLAND Collaboration], |arXiv:hep-ex/0205041| . 

[10] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90, 021802 (2003) 
|arXiv:hep-ex/021202"T] . 

[11] M. C. Gonzalez-Garcia, M. Maltoni and T. Schwetz, Phys. Rev. D 68, 053007 (2003) 
|arXiv:hep-ph/03062"26| . 

[12] The MINOS Collaboration, "Proposal for a Cosmic Ray Veto Shield for the MINOS 
Far Detector", Fermilab P-928, June 2002. 

[13] M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B 466, 415 (1999). 
|arXiv:hep-ex/9907037| . 

[14] Y. Declais et al, Nucl. Phys. B 434, 503 (1995). 

[15] F. Dydak et al, Phys. Lett. B 134, 281 (1984). 

[16] G. Barenboim, A. Dighe and S. Skadhauge, Phys. Rev. D 65, 053001 (2002) 
|arXiv:hep-ph/01060"02| . 

[17] S. Skadhauge, Nucl. Phys. B 639, 281 (2002) | [arXiv:hep-ph/0112189 . 

[18] A. Strumia, Phys. Lett. B 539, 91 (2002) |arXiv:hep-ph/0201134| . 

[19] Todd Haines, private communication. 

[20] M.S. Chanowitz, |arXiv:hep-ph/0304199| 

[21] T. Shutt [Borexino Collaboration], Nucl. Phys. Proc. Suppl. 110, 323 (2002). 

[22] The various solar neutrino scenarios are reviewed in K. Hagiwara et al. [Particle Data 
Group Collaboration], Phys. Rev. D 66, 010001 (2002). 



13 



P. C. de Holanda and A. Y. Smirnov, arXiv:hep-ph/0307266, 



14 



sub-GeV electrons 



u 



cos 



multi-GeV electrons 



cos 



sub-GeV muons 



COS 



multi-GeV muons 



i u 



cos 



Figure 3: SK zenith angle distributions normalized to no-oscillations expectations, for our 
first CPT violating example. Circles with error bars correspond to SK data. 
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Figure 4: SK zenith angle distributions normalized to no-oscillations expectations, for our 
second CPT violating example. Circles with error bars correspond to SK data. 
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